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Abstract 


Within this study a novel high fuel utilizing (High-uf) SOFC system is presented with special focus on the formation of nickel oxide, system 
efficiency and the required cell area at a fixed system performance of 1 MW. Within the High-uf SOFC cycle, a second SOFC stack is used to utilize 
a further part of the residual hydrogen of the first SOFC stack. This could be feasible by using an anode gas condenser, which is implemented 
between the first and the second stack. This reduces the water fraction of the anode gas and thereby the tendency of nickel oxide formation in case 
of a further fuel utilization. Thus, a higher total fuel utilization can be reached with the second SOFC stack. 

With the High-uf SOFC cycle, the system efficiency is increased by 7%-points compared to the simple atmospheric SOFC cycle. Furthermore, the 
average cell voltage and the fuel utilization are varied to carry out a first optimization of the stack’s power density. The results of this optimization 
have shown that the required cell area of the simple SOFC cycle can be slightly reduced by decreasing the fuel utilization, whereas the High-uf 
SOFC cycle shows an opposite effect. Here, the required cell area can be reduced at constant voltages by increasing the fuel utilization. Thus, 
higher system efficiencies could be reached with the High-uf SOFC cycle by using the same cell area as the simple SOFC cycle and at the same 
tendency of nickel oxide formation. 

A second condenser behind the second SOFC stack could be used to increase the carbon dioxide mass fraction up to 92%. This could be 


interesting for CO-sequestring applications as well. 
© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Fuel cell systems have the potential to reduce the fuel con- 
sumption, noise and gaseous emissions significantly [1]. The 
commercialization of fuel cell systems needs further develop- 
ments in materials, stack power density, stack durability and 
system design. These issues are strongly related with each other. 
The choice of electrochemical, thermodynamic and constructive 
parameters, which influence the system cost, have to be inves- 
tigated carefully for each specific application and each specific 
economic boundary condition. 

Previous studies have shown that solid oxide fuel cells 
(SOFC) coupled with gas turbines promise highest system effi- 
ciencies [1,2]. This is mainly caused by the comparably low 
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exergetic losses within the SOFC stack [3]. Siemens West- 
inghouse demonstrated the high durability of an atmospheric 
(system pressure ~ 1.3 bar) SOFC system with a performance 
of 100 kW. The efficiency of such atmospheric SOFC systems 
mainly depends on the average cell voltage, the electrochemical 
fuel utilization and the demand of excess air. The fuel utilization 
is commonly chosen with 85% at Ni-cermet anodes. A further 
amount in fuel utilization would result in a stronger formation 
of nickel oxide, which decreases the catalytic activity for the 
hydrogen oxidation. This effect is partly reversible but has to be 
avoided to increase availability and durability. 

Within this study a high fuel utilizing atmospheric SOFC 
system is presented with special focus on the formation of nickel 
oxide, system efficiency and the required cell area at a fixed 
system performance of 1 MW. This system is compared with 
a simple atmospheric SOFC system to estimate the benefits in 
efficiency and required cell area, which is strongly related to the 
installed stack cost. In this context, simulation models including 
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Fig. 1. Control volume. 


the SOFC, methane reformer and other system components have 
been developed within Matlab Simulink. 


2. Model 
2.1. Fuel processor mode 


For the simulation of the fuel processor, the thermodynamic 
equilibrium is assumed to determine the equilibrium composi- 
tions. The chemical species CH4, H2, H20, CO and CO%2 are 
taken into account. It is assumed that all gases are ideal fluids. 
The considered chemical reactions inside the reformer are the 
steam reforming 


CaHg +aH,0 < (£ + a) Hz + &CO 


and the shift reaction 
CO + H20 < H+ CO2 


The thermodynamic equilibrium of each reaction is given by 


K oy An A.G 
ES = = (=) exp (- : ) (1) 
p p RT 


where x is the molar fraction of reacting species i, p the abso- 
lute pressure, T the absolute temperature, A,G the free reaction 
enthalpy, A,H the reaction enthalpy and A,S is the reaction 
entropy. To avoid thermodynamically the carbon deposition (car- 
bon activity < 1), the steam/carbon ratio is chosen with 2.2. 


2.2. SOFC model 


The thermodynamic equilibrium of the steam reforming and 
the shift reaction is also used to determine the equilibrium com- 
positions inside the SOFC. At the anode, CH4, H2, H20, CO, 
CO? and at the cathode, O2, N2 are considered as species. The 
considered reactions inside the SOFC are: 

CHy + H20 < 3H) + CO 
CO + H20 <> CO2 + H2 
Hz + 0.502 <> H2O 


steam reforming : 

Shift Reation : 

Hydrogen Oxidation : 
The electrochemical oxidation of carbon monoxide is 


neglected. If the shift reaction fulfills the equilibrium condi- 
tion, the reversible work of the hydrogen oxidation and carbon 


monoxide oxidation will have the same value. Investigations of 
the shift reaction at Ni-cermet anodes have shown [4,5] that the 
conversions are always close to the equilibrium compositions. 

The conversion of the hydrogen and oxygen is determined by 
Faradays law 


I ; ; I 
AF’ ñm, =m = 5p (2) 
where the index “0” is used for the molar fraction at the entry 
of the control volume (Fig. 1). 
The electrochemical fuel utilization uf inside a SOFC without 
any leaks can be calculated by the cell current Z and the inlet gas 
composition: 


03,0 — NO, = 


I 
uf = 
F(8ncH,,0 + 27c0,0 + 2ñp,,0) 


(3) 


Polarization is acommon parameter in the analysis of fuel cell 
performance. To calculate the relationship between cell voltage 
Ucen and current density i, a constant area specific resistance R 
is assumed 


Uce = UNernst — i = UNernst — Ri (4) 


If the electrical performance of the cell Pe), the average cell 
voltage Ucey and the fuel utilization uf are requested parameters, 
the required cell area is calculated by the average current density 
i as follows: 


I Pal 


A=== > 
i Ucen i 


(5) 


The Nernst voltage is considered to be the driving force. It 
depends strongly on the partial pressures of reactants 


ArGH(T, P) 


UNemst = — OF 


PH20 


PH) 4/ PO 


Thus, the local current density changes along the flow direc- 
tion of the SOFC. This change is taken into account by the 
integral of the inverse overpotentials in the following equation: 


1 R uf2 1 
: f duf (7) 
u: 


i uf Juafı UnNernst(uty — Ucen 


1 
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Fig. 2. SOFC area map at design point. 


The results of Eq. (5) are shown in Fig. 2. The exemplary 
required cell area strongly depends on overpotential, fuel uti- 
lization and average cell voltage. The efficiency is defined 
as quotient of the total power and the enthalpy of reaction. 
The required cell area at a fixed performance decreases with 
a decreasing cell voltage. As a result the fuel consumption 
increases due to the lower efficiency of the fuel cell. This is 
an opposite effect which consists of a design optimum. 


2.3. Nickel oxide formation 


The formation of nickel oxide at the anode of a SOFC 
decreases the catalytic activity for the hydrogen oxidation. The 
formation of nickel oxide is partly reversible but has to be 
avoided to maximize the long-term stability at high fuel utiliza- 
tion. The oxidation of nickel could occur due to the reaction with 
H20 and CO. If the oxidation of hydrogen and carbon monoxide 
comply the equilibrium condition including the shift reaction, 
the oxygen activity can be calculated by one of these reactions. 

Different models were supposed for the electrochemical oxi- 
dation of hydrogen at Ni-YSZ (yttria stabilized zirconia) cermet 
anodes. The major differences of the models were found with 
regard to the location where the chemical and electrochemical 
reactions occur. De Boer suggested that interstitial hydrogen 
and hydroxyl are formed [7]. Bieberle as well as Mizusaki et al. 
suggested the presence of adsorbed oxygen on the Ni [6]. Mod- 
eling results of the anodic surface coverages of the species O, 
H, OH and H20 indicated that the fraction of adsorbed oxygen 
increases drastically with higher overpotentials [6], which could 
be an evidence of the formation of NiO at higher overpotentials. 

As long as the detailed mechanisms of the hydrogen oxida- 
tion and formation of nickel oxide at nickel cermet anodes are 
not yet clarified, the following assumptions are made for this 
work. At OCV conditions (open circuit voltage) the formation 
of NiO can be estimated by the equilibrium activity of the oxy- 
gen regarding the H2O—H) ratio of the gaseous phase. In case 
of operation and non-equilibrium condition at the anode, the 
potential of the anode (versus air) is displaced to less negative 
values, mainly caused by the activation overpotential. As a result 
of the following calculation, the concentration overpotential is 
lower than 10 mV. Thus, the concentration-dependent change in 


oxygen activity is considered to have a minor influence on the 
formation of NiO. The formation of NiO was observed at nega- 
tive anodic potentials between —850 and —650 mV (versus air) 
[7], which is the same range of oxygen activity governed by the 
equilibrium of Ni—-NiO system. This is why the anodic potential, 
calculated by activation and concentration versus air, is used in 
this work to estimate the formation of NiO. 

The formation of nickel oxide is estimated by comparing the 
anodic potential with the oxygen activity governed by equilib- 
rium of the Ni—NiO system. The equilibrium oxygen activity of 
Ni + 1/202 <= NiO is calculated by the following equation: 


log(a0,);_Nio = 8-96 — ——— (8) 


where po, ,Ni-NiO = 40,,Ni-Nio:| bar. To avoid the formation 
of nickel oxide at electrochemical equilibrium condition, it is 
assumed that the oxygen activity, given by the hydrogen oxida- 
tion (Eq. (9)), has to be lower than the oxygen activity of the 
nickel oxidation (Eq. (8)). The equilibrium of both equations 
(Eqs. (8) and (9)) were determined by measurements [8]. The 
oxygen activity of H? + 1/202 <> H20 is given by [8] 


26000 
Pi) -L 4594 (9) 
PH T 


log(ao, )H -H20 = 2log ( 


where py, is the partial pressure of the hydrogen and pp,o 
is the partial pressure of the water. The ratio of these partial 
pressures is strongly related to the fuel utilization of the SOFC. 
The results of this empirical approach is similar to those of the 
thermodynamic equilibrium. 

To compare the equilibrium of the Ni-NiO system with the 
overpotential of an anode, the oxygen activity of the nickel 
oxidation is converted into a electrical potential related to the 
potential of a reference air electrode (p= 1 bar and T= Tsofc). 
The anode potential versus air is given by 


AUa~n/Air = AUR, —H20/Air + fan 
= AUS, -m0 JAir + ÎAct,An + ÎDiff,An (10) 


where nan is the overpotential of the anode, taking the acti- 
vation (Act) and diffusion overpotential (Diff) of the hydrogen 
oxidation into account. The calculation of the activation and con- 
centration overpotential is described in the following Sections 
2.3.1 and 2.3.2 


a/ PO, Ait ) aD 
a/ PO2,Ni—NiO 


The equilibrium potential of the anode related to the reference 
air electrode can be determined by 


„/PO;,Air ) 
J Poz, H2—H20 


RT i re Pe (12) 
2F Pmo v p° 


Fig. 3 shows the tendency of the nickel oxidation (grey area), 
the potential equivalent of the hydrogen oxidation (dashed lines) 


r RT 
AUŅi-NiO/Air = oF In 


ij RT 
AUR, —H0/Air = oF In 


Q 
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Fig. 3. Anode potential (H2—H20) vs. air (1 bar). 


and the anode potential (full lines) as a function of the temper- 
ature and py, / ppo ratio. 

The formation of nickel oxide thermodynamically occurs at 
oxygen pressures higher than 5.9 x 107!! bar (1000°C). This 
partial pressure is equal to an anode potential of —603 mV ver- 
sus air (1 bar, 1000 °C). Guindet et al. observed an influence on 
nickel oxide formation at an anode potential, more positive than 
—800 mV (950°C) [9]. Therefore, the anode potential is chosen 
lower than —800 mV in this study. This value is reached at a 
PH / pmo ratio of 0.15 (fuel utilization = 85%) and at a temper- 
ature of 1000 °C. This is particularly important at the gas outlet 
of the anode where high fuel utilizations occur. The influence 
of the activation and diffusion overpotential on the anode poten- 
tial is comparably small in this range of temperature and fuel 
utilization. The activation overpotential increases with lower 
temperatures and the diffusion overpotential increases mainly 
with lower py, / pmo ratios. 


2.3.1. Local activation overpotential 

The activation of a chemical reaction involves energy barriers 
which must be overcome by the reacting species. These energy 
barriers are called the activation energy and result in charge 
transfer polarization. The activation polarization of an electro- 
chemical reaction is usually expressed by the Butler—Volmer 
equation: 


Sele lex âz F Fact ex (1 — &)zF Fact (13) 
as [P| RT P RT 


where @ is the transfer coefficient and ip is the exchange cur- 
rent density. The transfer coefficient is approximately 0.5 for 
fuel cells. The exchange current density is the cathodic and 
anodic electrode reaction rate at the equilibrium potential. A 
high exchange current density results in a good performance. 
With â = 0.5 the activation polarization can be calculated as 


follows: 
RT i i 
arsinh (=) (14) 


AAct = ——— 
i 0.5zF io 


For the anodic exchange current density of the hydrogen 
oxidation at Ni-YSZ anodes, Mogensen [10] has proposed the 


following equation: 


. H H20 E Act, A 
io,an = YAn (=) (22) exp (-=44") (15) 
PO,An Po,An RT 


where Fact,an is the activation energy, yan an empirical pre- 
exponential factor, py, the partial pressure of hydrogen, py,o 
the partial pressure of water and po,An is the reference pressure 
at the anode. 


2.3.2. Local concentration overpotential 

The concentration overpotential at electrodes is a fictive elec- 
trical potential drop caused by the concentration difference in the 
gaseous phase. Diffusion through porous materials is typically 
described by the binary molecular diffusion and the Knud- 
sen diffusion. In this study, both types of diffusion are taken 
into account to estimate the concentration overpotential. It is 
assumed that the concentration gradient along the pore is much 
higher than the concentration gradient within the surface bound- 
ary. Thus, the concentration gradient inside the surface boundary 
is neglected. 

For the ideal gas model the anodic concentration overpoten- 
tial, which is proportional to the irreversible entropy production, 
can be written as 


RT ty PE 
ADitf,An = -y In Ph, RO (16) 
Z PH, Pho 


where p° is the partial pressure in the gas bulk and p’ is the partial 
pressure at the end of the pores. The concentration profiles inside 
the pores can be estimated by using the Fick’s law. The molar 
fraction x can be expressed as 


RTiô 

Paleo (17) 
: 2 pP2FDeff An 
RTiôAn 


(18) 


* — 0 
ia 7 *H20 a p2FDeff, An 
for the equimolar diffusion of hydrogen and water. dan is here 
the thickness of the anode and Deff is the effective diffusion 
coefficient. The effective diffusion coefficient takes the Knudsen 
diffusion as well as the binary diffusion at the anode into account. 
The binary diffusion coefficient between gas species A and 
B in free space is given by the Chapman—Enskog theory 
1 1 ) ua g3 


Dap = 1.8583 x =( ae 


19 
Ma | Mp (19) 


P TAB ND 
where o is the collision diameter (Lennard-Jones length) and Qp 
is the collision integral based on the Lennard-Jones potential. 
The collision integral is calculated by the Lennard-Jones energy 
€. The Lennard-Jones energy is a function of the critical gas 
properties and the Boltzmann constant k [11]. 

The Knudsen diffusion coefficient for a species A is defined 
as 

2Rporda 


Dxkn,A = a (20) 


where Rpor is the pore radius and J is the average velocity of the 
gas molecule. Both ordinary diffusion and Knudsen diffusion 
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can be combined within the effective diffusion coefficient which 
can be written as 


i t / 1 1 
= + (21) 
Deff Pora 9 \ Dap Dkna 


The tortuosity Tp and porosity @ are used to account the geom- 
etry of the pores and relation between solid space and pore space 
of the electrode. 

If the presence of other species except hydrogen and water 
inside the anodic pores is neglected, the anodic effective diffu- 
sion coefficient can be expressed for the ideal gas-mixture as 
[12]: 


PH20 PH 
Detf,An = ( Z ) Deff, Por, H> + (=) Deff,Por, mo (22) 
Pges Pges 


(Table 1) lists the assumptions which are made for the previ- 
ous estimation of the anodic overpotentials. 


3. System configuration 
3.1. Simple SOFC cycle 


The configuration of the simple SOFC cycle with a perfor- 
mance of 1000 kW is shown in Fig. 4. The simulation is based 
on the solution of the energy and the mass balances calculated 
within Matlab Simulink. The model of each component consists 
of a mode for design studies and a off-design mode which takes 
the part load characteristics of the fuel cell, the turbo machinery, 
the pumps and the heat exchangers into account. 

The thermodynamic equilibrium Eq. (1) is used to solve the 
mass balance of the reformer for a predicted outlet temperature. 
In the reformer the hydrocarbons are converted to H2, H20, 


Quel (Hu) = 1882.5 KW 
Pel= 1000 kW 


miz 53.12 % 


Water Methane Air 


Table 1 

Anode parameters 

Thickness anode (Ni-YSZ) (m) 50 x 1076 
Pre-exponential coefficient for ig (A m~?) 18.9 x 10° 
Activation energy anode (kJ mol7!) 120 

Pore radius (m) 0.5 x 1076 
Tortuosity 6 

Porosity (%) 30 


CO, CO, and residual CH4 by combining the steam reforming 
and the shift reaction. For the SOFC stack a constant outlet gas 
temperature of 1000°C is chosen. To fulfil the energy balance 
of the SOFC stack, the inlet temperatures of the SOFC stack is 
calculated iteratively for a requested total excess air ratio. That 
involves a partial preheating of the reactants inside stack. 

The air is preheated in the heat exchanger (2b) and in a stack 
integrated heat exchanger to reduce the excess air demand. In the 
combustor (6a) the methane, hydrogen and carbon-monoxide of 
the anode gas is oxidized to use the residual heat. The evaporator 
(1b) is implemented in the exhaust gas stream. 

With the listed assumptions (Table 2) the simple SOFC cycle 
reaches an electrical efficiency of 53% (LHV) at a system 
pressure of 1.27 bar, an average cell voltage of 0.7 V and an 
electrochemical fuel utilization of 85%. 


3.2. High-uf SOFC cycle 


The configuration of the high fuel utilizing (High-uf) SOFC 
cycle with a performance of 1000 kW is shown in Fig. 5. This 
system consists of an additional SOFC stack which is connected 
in series with the first SOFC stack. The second SOFC stack is 
supplied by the cathode gas of the first SOFC stack after it passed 


Fig. 4. Simple atmospheric SOFC cycle. 
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Table 2 

Initial parameters of simple SOFC cycle 

SOFC average cell voltage (V) 0.7 

Fuel utilization (%) 85 

SOFC temperature gas outlet (K) 1273.15 

System pressure (bar) 1.27 

Air ambient conditions (dry air) 288 K, 1.011 bar 
Total excess air ratio 1.44 
Steam/carbon ratio 2.2 

Isentropic efficiency of blower (%) 87 

Mechanical efficiency of blower (%) 95 

Pressure drop SOFC, HEX, combustor (bar) 0.04; SOFC: 0.01 
Thermal losses SOFC, HEX, combustor (%) 0 

Inverter efficiency (%) 95 

Motor efficiency (%) 96 


the air preheater (2b). The anode gas of the first SOFC stack is led 
through the regenerative heat exchanger (1d) and the condenser 
(le) before it enters the second SOFC stack. The fuel utlization 
of the first SOFC stack is chosen with 85%. This is equal to 
the fuel utilization of the simple SOFC cycle. The second SOFC 
stack is used to increase the fuel utilization up to 97% at an aver- 
age cell voltage of 0.7 V. It is assumed that the maximum fuel 
utilization is limited by the py, / pmo ratio of partial pressures 
at the outlet of the anode (Section 2.3). The first stack reaches a 
PH / pmo ratio of 0.12 and a local Nernst voltage of 0.77 V at 
the anode outlet (Fig. 6). Within the condenser (le) the partial 
pressure of the water is reduced. Thus increases the py, /PH,o 
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Fig. 6. Nernst voltage vs. fuel utilizations. 


ratio up to 0.46 at the condenser outlet (le). If the equilibrium 
condition of the shift reaction is taken into account, the reaction 
will generate water and carbon monoxide by consuming hydro- 
gen and carbon dioxide at the inlet of the second SOFC stack. If 
the shift reaction complies their equilibrium, the py, / pH,o ratio 
will decrease again and a Nernst voltage of 0.94 V at the inlet of 
the second SOFC stack (Fig. 6) is reached. With this py, /PH,0 
ratio it is feasible to carry out a further fuel utilization without an 
amount in nickel oxide formation. The fuel utilization of the first 
and the second SOFC stack increases to 97% in sum when the 
PH>/ pmo ratio (0.12) reaches the same value as the first SOFC 
stack. Another advantage of this configuration is the amount in 
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Fig. 5. High-uf atmospheric SOFC cycle. 
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power density which is due to the higher Nernst voltage at the 
second SOFC stack inlet. Furthermore, the second condenser 
(1f) increases the carbon dioxide mass fraction to 92%, which 
could be interesting for CO2-sequestring applications. High car- 
bon dioxide mass fractions reduce the energy consumption of 
the CO storage. 

The High-uf SOFC cycle reaches an electrical efficiency of 
60% (LHV) at a system pressure of 1.27 bar, an average cell 
voltage of 0.7 V and a total electrochemical fuel utilization of 
97%. 

As mentioned before, the steam/carbon ratio is chosen above 
2 to avoid carbon deposition. This ratio is reduced within the 
condenser le. Even if the methane is approximately complete 
converted in the reformer and the first SOFC stack and no carbon 
deposition is expected from that specie, the Boudouard reaction 
could result in carbon deposition, especially at low tempera- 
tures. This issue is similar to the effort, supplying a PEFC by 
the anode gas of a SOFC. For this particular configuration, a 
shift reactor, located between the SOFC and the PEFC, is used 
to convert the residual carbon monoxide and water into carbon 
dioxide and hydrogen. A side effect of the shift reaction is, that 
the carbon activity due to the Boudouard reaction is reduced in 
that case. This could also become necessary for the proposed 
“High-uf SOFC cycle”, if a carbon deposition occurs in the 
condenser le. 


3.3. Required cell area 


To compare these two SOFC configurations, the required cell 
area (Eq. (5)) is estimated at a constant average cell resistance 
of 0.92 Qcm?. The reference case of optimization is chosen at 
Ucen = 0.7 V, uf =85% of the simple SOFC cycle. The system 
efficiency is defined as 


Pat,sys 


— (23) 
ACH,,sys LHV 


Nel,sys = 
where Pelsys is the system performance, ncHy,sys the molar 
methane consumption of the system and LHV is the molar lower 
heating value of the methane. It can be shown that the total elec- 
trical current (Eq. (5)) of parallel cells, is also influenced by the 
requested system’s performance and the system’s efficiency 


Pel,sys 


I = uf F8ncu, = uf F8 ——_-#!__ 
EMESIS Nel, sys LHV 


(24) 


Eqs. (5) and (24) are used to calculate the change in required 
cell Asorc/Asorc.o at a constant average cell resistance and a 
fixed system performance of 1 MW. 


Asorc _ Nel,sys,0 uf io 
ASOFC,0 Nelsys Ufo i 


(25) 


Fig. 7 shows the influence of the fuel utilization and the aver- 
age cell voltage on the change in the total required cell area 
of the simple SOFC cycle (doted lines) and the High-uf SOFC 
cycle (full lines). The fuel utilization of the simple SOFC cycle 
is varied from 65% to 85%. For the High-uf SOFC cycle the 
fuel utilization of the first SOFC stack is kept constant (85%). 
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Fig. 7. Influence of U and uf on the required cell area. 


The total electrochemical fuel utilization is varied from 85% to 
97% with the second stack. Fig. 7 shows that the same system 
efficiency can be reached with different combinations of cell 
voltage and fuel. Depending on this combination (Ucey, uf) the 
thermodynamic parameters of the process and the required cell 
area are changing. 

The average cell voltage of the SOFC is varied from 0.6 to 
0.75 V. The required cell area of the simple SOFC cycle can be 
slightly reduced by decreasing the fuel utilization because the 
average power density of the stack increases with lower fuel 
utilizations (Fig. 8). The required cell area of the High-uf SOFC 
cycle shows an opposite effect. Here, the required cell area can 
be reduced at constant voltages by increasing the fuel utilization. 
This is mainly caused by the amount in system efficiency and 
average power density of the second stack (Fig. 8). The high 
hydrogen fraction at the inlet of the second stack increases the 
power density significantly. In case of the High-uf SOFC cycle 
the influence of the system efficiency on the area reduction is 
stronger than the influence of the fuel utilization on the area 
increment. 

The average power density of the fuel cell (Fig. 8) is calcu- 
lated by the stack’s total power and the total required cell area 
at a constant average cell resistance of 0.92 Q cm”. 

With the presented High-uf SOFC cycle the system efficiency 
can be increased significantly by using the same cell area as the 
simple SOFC cycle. This could reduce the installed stack cost 
at a fixed system performance. 
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Fig. 8. Influence of U and uf on the average power density. 
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4. Conclusion 


Within this study a high fuel utilizing atmospheric SOFC 
system is presented with special focus on the formation of nickel 
oxide, system efficiency and the required cell area at a fixed 
system performance of 1 MW. In the High-uf SOFC cycle, a 
second SOFC stack is used to utilize a further part of the residual 
hydrogen of the first SOFC stack. This is feasible because the 
tendency of nickel oxide formation at the anode decreases with 
a lower water fraction. To reduce the water fraction of the anode 
gas, acondenser is implemented between the first and the second 
stack. With this configuration, a system efficiency of 60% is 
reached at a total (first and second stack) electrochemical fuel 
utilization of 97% and an average cell voltage of 0.7 V. This is 
7%-points higher than the efficiency of the simple atmospheric 
SOFC cycle. 

To compare these two configurations, the required cell area 
is calculated. The required area gives an idea of the expected 
change in stack cost and stack weight, respectively. Further- 
more, the average cell voltage and the fuel utilization are varied 
to carry out a first optimization of the stack design. The starting 
point of this optimization is chosen at an average cell voltage 
of 0.7 V and a fuel utilization of 85%. For the High-uf SOFC 
cycle the fuel utilization of the first SOFC stack is kept constant 
(85%). The total electrochemical fuel utilization is varied from 
85% to 97% with the second stack. The results of this optimiza- 
tion have shown that the required cell area of the simple SOFC 
cycle can be slightly reduced by decreasing the fuel utilization, 
whereas the High-uf SOFC cycle shows an opposite effect. Here, 
the required cell area can be reduced at constant voltages by 
increasing the fuel utilization. Commonly, the required cell area 
increases with an increasing average cell voltage and an increas- 
ing fuel utilization. But in case of the High-uf SOFC cycle the 
influence of the efficiency on the area reduction is stronger than 
the influence of the fuel utilization on the area increment. With 
the presented High-uf SOFC cycle the system efficiency can be 
increased significantly by using the same cell area as the simple 
SOFC cycle. 


A second condenser behind the second SOFC stack can be 
used to increase the carbon dioxide mass fraction up to 92%. 
This could be interesting for CO2-sequestring applications. 

Compared to the simple SOFC cycle, it is conceivable that 
the High-uf SOFC reduces the system cost as long as the cost 
of the additional heat exchanger and condenser is lower than 
the savings in stack cost. Furthermore, the increment in system 
efficiency allows higher system cost. But this depends strongly 
on the application. 

Future work will involve more detailed considerations of the 
high fuel utilizing stack with regard to the local temperature and 
current distribution of the cell. 
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